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ABSTRACT 


Frequency  domain  methods  are  employed  for  the  design  and 
analysis  of  anti-roll  tank  systems.  A  waveslope  spectrum  of 
the  ocean  is  developed  based  upon  the  trochoidal  model  of 
ocean  wares  and  the  Pierson-Moskowitz  ocean  wave  power 
spectra.  The  response  of  a  typical  ship  to  the  waveslope 
spectrum  is  studied  with  the  aid  of  a  numerical  example. 

FORTRAN  subroutines  for  the  design  of  passive  and 
nump-activated  O-tube  anti-roll  tanks  are  developed,  and  the 
effects  of  a  tank  installation  on  the  ship  parameters  is 
observed.  Many  controller  design  methods  are  investigated, 
and  final  design  is  accomplished  by  exhaustive  search  of 
closed-loop  frequency  responses.  Frequency  responses  of 
passive  and  activated  tank  stabilized  ships  are  compared 
with  that  of  an  unmodified  ship  to  illustrate  the 
effectiveness  of  design.  The  effects  of  vessel  displacement 
variations  due  to  consumption  of  fuel  and  stores  on  tank 
effectiveness  and  controller  design  are  discussed. 
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I. 


FORMULATION  OF  THE  PROBLEM 


A  ship  on  the  ocean  is  in  constant  motion.  Ideally,  the 
only  motions  would  be  translation  in  the  forward  direction 
and  rotation  about  the  vertical  axis,  associated  with 
propelling  and  maneuvering  the  ship.  Realistically,  there 
are  perturbation  motions  in  addition  to  the  desired  motions 
in  the  form  of  translations  along  and  rotations  about  the 
longitudinal,  transverse,  and  vertical  axes  of  the  ship. 
The  purpose  of  this  thesis  was  to  study  the  rotation  about 
the  longitudinal  axis,  known  as  roll,  and  arrive  at  a  scheme 
for  controlling  it.  A  procedure  for  designing  the 

installation  of  a  U-tube  anti-roll  tank  system  was 

consolidated  from  several  sources  ["chadwick,  1 950~ )  ["chadwick 
and  blotter,  1 

1973^  [ziegler ,  1 

subroutine.  Many  controller  design  methods  were 

investigated,  and  the  final  design  was  accomplished  by 

exhaustive  search  of  closed-loop  frequency  responses. 

A.  HAVSSLDPE  SPECTRUM  OF  THE  OCEAN. 

Before  a  meaningful  analysis  of  ship  motions  in  water 
could  be  conducted,  the  motions  of  the  water  itself  had  to 
be  understood.  Because  a  ship  tends  to  ride  normal  to  the 
surface  of  the  water,  the  most  significant  physical  property 
of  waves  with  respect  to  ship  motions  is  the  waveslope, 

At  the  crest  and  in  the  trough  of  a  wave,  the  waveslope  is 
zero,  while  at  some  point  nearly  midway  between  those 
extremes,  the  waveslope  is  maximum. 

Several  theories  of  wave  motion  have  been  advanced, 
including  the  simple  harmonic  wave;  the  Stoke's  Nave,  which 


953^j  [^Blagoveshchensky,  1962j|  Vachanaratana , 
9751  and  presented  in  the  form  of  a  FORTRAN 
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is  a  harmonic  wave  with  a  second  harmonic  added;  and  the 
trochoidal  wave.  Geometrically,  a  trochoidal  curve  is 
constructed  by  rolling  a  circle  of  radius  R  under  a 
horizontal  line,  and  tracing  the  locus  of  a  point  a  radial 
distance  c<?.  from  the  center  of  the  circle.  The  trochoidal 
curve  may  be  steep  or  shallow  depending  on  the  ratio  of  r:R, 
and  corresponds  very  closely  to  actual  observed  ocean  waves, 
where  2  ir  R  =  L  =  wavelength,  and  2r  =  h  =  trough- to-crest 
wave  height. 

Ocean  waves  are  primarily  caused  by  the  transfer  of 
kinetic  energy  from  the  wind,  which  is  the  result  of  the 
transfer  of  thermal  energy  in  weather  phenomena.  As  the 
wind  blows  over  a  flat  surface  of  water,  the  water  particles 
at  the  air-sea  interface  are  excited  by  frictional  forces, 
and  tend  to  move  downwind.  The  freedom  of  motion  of  these 
particles  is  limited  by  the  fluidic  properties  of  water,  and 
thus  they  tend  to  accumulate  in  successive  ridges,  or 
ripples,  on  the  sea  surface.  Once  these  ridges  are  formed, 
further  wave  buildup  continues  by  a  different  mechanism. 
On  the  upwind  side  of  a  ridge,  a  small  local  high  pressure 
region  develops  which  exerts  a  downward  force  on  the  water 
surface,  forming  the  wave  trough.  On  the  downwind  side  of  a 
ridge,  a  corresponding  small  local  low  pressure  region 
develops  which  draws  the  water  surface  upward,  forming  the 
wave  crest.  The  wave  continues  to  grow  as  long  as  the 
wind  continues  to  supply  energy  in  excess  of  that  consumed 
by  the  internal,  friction  of  the  water  in  the  wave,  which  it 
does  until  the  wave  velocity  equals  the  wind  velocity. 

Three  key  factors  govern  the  development  of  wave 
systems.  They  are  the  initial  difference  between  the 
velocities  of  the  wind  and  any  existing  wave  systems;  the 
length  of  time  (duration)  that  the  wind  continues  at  a  given 
velocity;  and  the  fetch,  or  distance  of  relatively  uniform 
depth  water  over  which  the  wave  will  run  before  being 
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altered  by  a  change  of  depth  or  destroyed  by  encountering  a 
land  mass.  (The  passage  of  a  vessel  through  an  ocean  wave 
system  does  little  to  alter  the  growth  or  shape  of  the 
system . ) 

The  rate  of  increasing  wave  height  is  a  function  of  the 
initial  difference  between  the  velocities  of  the  wind  and 
the  existing  wave  system.  A  sudden  change  of  wind  velocity, 
such  as  that  associated  with  the  passage  of  a  front,  can 
cause  the  waves  to  increase  at  a  rate  of  one  or  two  feet  per 
minute.  There  is  a  physical  limit  to  the  increase  in  wave 
height,  and  it  is  defined  in  terms  of  the  wave  height  to 
wave  length  ratio,  h:L=1:7.  Near  this  ratio  limit,  excess 
wind  enecgy  blows  off  the  tops  of  the  waves,  forming 
turbulence  known  as  "white  caps." 

The  rate  of  increasing  wave  length  generally  lags  behind 
the  rate  of  increasing  wave  height.  Thus  in  regions  of 
short  fetch,  a  sudden  change  of  wind  velocity  may  never 
produce  waves  with  long  wavelengths  due  to  the  limiting 
ratio.  In  contrast,  on  the  open  ocean  where  there  is  a  long 
fetch,  a  change  of  wind  velocity  can  produce  long  waves  of 
relatively  low  height.  For  waves  of  length  L>100  feet,  there 
is  a  diminishing  likelihood  of  ever  approaching  the  maximum 
value  of  the  ratio.  In  practice,  for  wavelengths  between 
500  and  500  feet,  the  ratio  rarely  exceeds  1:20;  for  the 
longer  ocean  waves  (of  length  over  1000  feet),  the  ratio 
drops  to  1:50. 

The  above  discussion  of  wave  development  only  hints  at 
the  complexity  of  attempting  to  define  the  waveslope 
spectrum  of  the  ocean.  Conferences  have  been  held,  and  many 
books  have  been  written  on  the  subject  of  ocean  spectra. 
(^Conference,  Easton,  Maryland,  196lJ  j_Inoue,  1  966J 
(^Comstock,  1  96 7 j  ^Gillmer,  196 9~ |  fziegler,  1975j  Many  ocean 
spectra  plot  relative  wave  energy  as  a  function  of  frequency 
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with  "sea  state"  as  a  parameter.  Others  define  wind 
velocity  and  duration  as  parameters.  Still  others  define 
"significant  wave  height"  as  a  parameter,  with  many 
definitions  of  "significant  wave  height"  itself.  Spectra 
have  been  constructed  for  given  locations  citing  mean  wave 
height,  or  mean  squared  wave  height,  or  RUS  wave  height  as 
parameters  . 

Two  other  factors  which  affect  the  effective  ocean 
spectra  are  vessel  speed  and  angle  of  encounter.  To 
illustrate  these  effects  on  the  ocean  spectrum,  consider  a 
ship  proceeding  with  velocity  identical  to  that  of  a 
single-frequency  wave  train.  If  the  ship  heads  directly 
into  the  approaching  wave  train,  the  apparent  frequency  of 
the  wave  doubles,  but  the  roll  coupling  with  the  ship  goes 
to  zero.  When  the  ship  travels  parallel  to  the  wave  train, 
which  is  to  say  the  ship  is  in  a  beam  sea,  the  apparent 
frequency  of  the  wave  is  the  same  as  the  actual  frequency  of 
the  wave,  and  the  roll  coupling  is  maximum.  when  the  ship 
heads  in  the  same  direction  as  the  wave  train  ("following 
seas"  condition),  the  apparent  frequency  of  the  wave  drops 
to  zero,  as  does  the  roll  coupling.  If  B  denotes  the  ship's 

O 

angle  of  encounter  with  the  wave  train,  where  B=0  for  a 

.  O  O 

following  sea,  B=90  for  a  beam  sea,  and  B=180  for  a  head 
sea,  then  the  apparent  waveslope  is  related  to  the  actual 
waveslope,  and  the  apparent  frequency  of  the  sea  is  related 
to  the  actual  frequency  by  the  expressions  Pziegler,  197 5~1 : 


*<B)  = 

4>  sin  B 

max 

(1) 

w  (B)  = 

w  (l.-u)  u  /g) 

sea  sea  ship 

cos  B 

(2) 

ship 

forward  velocity  of  the 

ship 

g  =  acceleration  due  to  gravity 
Since  the  ultimate  reason  for  this  study  was  to  aid  in 
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the  desigi  of  a  general-purpose  anti-roll  stabilizer,  the 
ocean  spectrum  chosen  was  a  wave  height  power  density 
spectrum,  with  sea  state  and  significant  wave  height  (1/3 
highest  waves)  as  parameters.  The  waveslope  spectrum  was 
then  computed  on  the  basis  of  the  trochoidal  model  of  the 
ocean,  where 


r 

1/f  = 

2  r  /  co 

seconds 

(3) 

L 

5.118 

2 

T 

feet 

(4) 

P  = 

max 

(h/L) 

X  57.3 

degrees 

(first  approximationj 

B.  SHIP  ROLLING  IN  WAVES. 

A  ship  on  the  ocean  moves,  in  general,  with  six  degrees 
of  freedom  relative  to  the  coordinate  axis  system.  By 
convention,  the  positive  x-direction  is  forward  along  the 
longitudinal  centerline  of  the  ship,  the  positive 
y-direction  is  horizontally  outward  on  the  starboard  beam, 
and  the  positive  z-direction  is  vertically  downward.  The 
center  of  the  coordinate  axis  system  is  located  at  the 
vessel  center  of  gravity.  Rotation  about  the  x-axis  is 
known  as  roll,  with  positive  roll  angle,  <t> ,  defined  as  an 
inclination  to  starboard.  Roll  can  be  considered  separately 
from  the  other  degrees  of  freedom,  which  reduces  the  study 
of  ship  motions  to  a  two-dimensional  problem  with  one  degree 
of  freedom. 


1.  Definitions  and  Concepts  of  Naval  Architecture. 


The  following  definitions  and  concepts  of  naval 
architecture  ^Comstock,  1967^j  Qiillmer,  1969^  are  given  for 
the  benefit  of  the  reader,  and  are  illustrated  in  Figure 

1-1 . 

a.  Displacement. 

The  displacement  of  a  vessel,  d ,  is  defined  as 
the  weight  of  the  volume  of  water  displaced  by  the  vessel. 
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G  =  Center  of  Gravity 
B  =  Center  of  Buoyancy 
M  =  Metacenter 
GM  =  Metacenter  Height 
GZ  =  Righting  Arm 

4)  =  Ship  Roll  Angle  Relative  to  Earth  Reference 

WL  =  Water  Line 


Figure  1-1.  Definitions  and  Concepts  of  Naval  Architecture 


Dua  to  the  Principle  of  Archimedes,  the  displacement  equals 
the  weight  of  the  vessel  itself.  Displacement  is  always 
given  in  long  tons,  where  one  long  ton  equals  2240  U.S. 
pounds.  The  weight  of  the  water  displaced  equals  the  weight 
of  the  vessel,  independent  of  the  density  of  the  water.  The 
volume  of  water  displaced,  however,  is  a  function  of  the 
density  of  the  water.  For  purposes  of  most  calculations, 
the  density  of  fresh  water  is  given  as  36  cubic  feet  per 
ton,  whereas  the  density  of  salt  water  is  given  as  35  cubic 
feet  per  ton . 

b.  Center  of  Gravity. 

The  center  of  gravity,  G,  is  defined  as  the 
point  at  which  the  summation  of  the  moments  of  all  the 
weights  in  a  vessel  with  reference  to  any  axis  passing 
through  that  point  equals  zero.  The  center  of  gravity  of 
most  warships  of  usual  form  is  generally  located  near  the 
designer's  waterline,  in  or  near  the  plane  of  the  midship 
section.  For  purposes  of  buoyancy  and  stability  analysis, 
the  entire  weight  of  a  vessel  is  considered  to  be 
concentrated  at,  and  the  force  of  the  weight  acts  vertically 
downward  through,  the  center  of  gravity. 

c.  Center  of  Buoyancy. 

The  center  of  buoyancy,  B,  is  defined  as  the 
geometric  canter  of  the  underwater  form  of  the  vessel.  The 
force  of  buoyancy,  equal  in  magnitude  to  the  force  of  the 
displacement  of  the  vassal,  acts  vertically  upward  through 
the  center  of  buoyancy.  When  a  vessel  is  in  static 
equilibrium,  the  center  of  buoyancy  is  located  directly 
under  the  center  of  gravity,  such  that  the  force  vectors  of 
displacement  and  buoyancy  are  equal  and  opposite.  Due  to 
this  definition  of  the  center  of  buoyancy,  as  the  condition 
of  a  vessel  varies  from  that  of  static  equilibrium,  the 
location  of  B  varies  also. 
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d.  Metacenter. 

As  a  ship  inclines,  the  shape  of  the  underwater 
form  changes,  and  the  center  of  buoyancy  migrates  to  remain 
at  the  geometric  center.  The  migrations  form  a  locus  which 
is  dependent  upon  and  directly  related  to  the  shape  of  the 
ship's  hull.  This  locus  is  generally  elliptic  and  as  such 
cannot  have  a  fixed  center,  but  for  an  infinitesimal 
increment  of  the  locus  the  instantaneous  center,  called  a 
metacenter,  can  be  located.  The  locus  of  metacenters  is 
geometrically  an  evolute,  and  is  called  the  metacentric. 
For  small  angles  of  inclination,  up  to  10°  for  most  ships, 
the  metacenters  are  co-located  at  a  point  known  as  the 
Metacentec,  M. 

e.  Metacentric  Height. 

When  a  ship  is  upright,  that  is,  when  the  angle 
of  inclination  <h  =  0°  ,  the  centers  M,  G,  and  B  lie  in  a 
straight  vertical  line  over  the  keel.  The  distance  between 
G  and  M  is  known  as  the  metacentric  height,  Gti.  If  M  is 
above  G,  GM  is  defined  as  being  positive. 

f.  Righting  Arm. 

When  a  ship  is  inclined,  the  center  of  buoyancy 
is  no  longer  on  the  original  vertical  centerline.  A  point  Z 
is  described  by  the  intersection  of  a  line  drawn  vertically 
upward  from  B  with  a  line  drawn  horizontally  away  from  G. 
The  distance  GZ  is  known  as  the  righting  arm,  and  for  small 
angles  of  inclination  <t>  , 

GZ  =  GM  Sin  <t>  =  GM  <t>  (6) 

g.  Righting  moment. 

The  righting  moment,  Mr  ,  is  the  product  of  the 
force  (displacement)  and  the  righting  arm: 

Mr  =  A  GZ  =  A  GM  <fi  (7) 

h.  Metacentric  Height  and  Stability. 

The  metacentric  height  is  a  measure  of  the 
statical  stability  of  a  vessel.  A  large  positive  GM 
corresponds  with  a  pair  of  high-f reguency  complex  conjugate 


15 


poles  in  tbe  left  half  of  the  s-plane,  far  removed  from  the 
j ^  -axis.  A  small  positive  GM  corresponds  with  a  pair  of 
low-frequency  complex  conjugate  poles  in  the  left  half  of 
the  s-plane,  close  to  the  jw  -axis.  If  GM  =  0,  then  there 
are  two  poles  at  tne  origin  of  the  s-plane.  A  negative  GM 
puts  poles  in  the  right  half  of  the  s-plane. 

2 .  Equations  of  Motion. 


a.  Natural  Period  of  Roll. 

In  still  water,  a  ship  can  be  made  to  roll  by 
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righting  moment  produces  rotation  of  the  ship  towards  the 
upright  position.  The  potential  energy  is  converted  to 
kinetic  energy  such  that  when  the  ship  is  in  the  upright 
position,  assuming  unresisted  roll,  all  the  potential  energy 
has  been  converted  to  kinetic  energy.  As  a  result,  the  ship 
will  continue  to  rotate  past  the  upright  position  until  all 
kinetic  energy  has  been  reconverted  to  potential  energy,  at 
the  opposite  limit  of  roll.  Assuming  no  loss  due  to 
friction,  the  ship  would  continue  to  oscillate,  or  roll, 
indefinitely  with  constant  amplitude.  The  period  of  these 
oscillations  can  be  found  from  simple  harmonic  motion 
relationships  to  be  [^Comstock,  1967^j 
T  =  (2 vk)  4/gGM 

=  (t .  103k)  /v/gF 

*  (0.44b)  /  sj  GM  (8) 

whe  re 


2 

g  =  acceleration  due  to  gravity,  ft/sec 
k  =  radius  of  gyration  of  mass  of  ship,  feet 
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b  -  extrema  beam  of  ship,  feet 


b.  Roll  by  Wave  Action. 

The  cause  of  ship  rolling  in  a  seaway  is 
primarily  the  unbalanced  moments  which  result  from  the 
center  of  buoyancy  shifting  in  response  to  the  change  of  the 
waterplane,  and  thus  the  shape  of  the  underwater  portion  of 
the  hull,  as  waves  pass  under  the  ship.  A  ship  will  usually 
attempt  to  ride  normal  to  the  surface  of  a  wave.  A  ship 
also  moves  in  a  circular  orbit  as  do  the  particles  of  water 
in  a  wave.  The  orbits  induce  centrifugal  forces,  which  are 
opposed  by  the  dynamic  forces  in  water. 

Rolling  in  waves  is  a  forced  oscillation,  and  the  ship 
will  tend  to  oscillate  at  or  near  the  frequency  of  the 
waves.  If  the  waves  are  not  uniform,  which  is  usually  the 
case,  a  ship  will  attempt  to  oscillate  with  its  natural 
frequency.  The  overall  roll,  therefore,  is  a  combination  of 
the  ship’s  natural  and  the  wave's  periods,  with  the  wave 
influence  being  generally  more  predominant. 

c.  Equations  of  Motion. 

A  ship  can  be  modeled  as  a  second-order, 
spring-mass-damper  system.  j^Chadvick,  1950J  ^Comstock,  1967J 
The  spring  constant,  or  static  righting  coefficient,  was 
discussed  earlier  in  this  paper  in  the  definition  of  the 
righting  moment.  The  static  righting  coefficient  is 

K  =  AGS  (9) 
s 

The  inertia  term,  or  roll  inertia  coefficient,  is  the 
product  of  the  mass  of  the  ship  and  the  radius  of  gyration 
of  the  mass  of  the  ship  squared: 


2 

J  =  (  A/g)  k  (10) 

s 

The  ndius  of  gyration  is  found  by  timing  the  natual 
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period  of  roll  in  protected  (still)  water,  and  massaging  the 
equation  for  roll  period  as  follows: 

k  =  (T+N/gTfT  )  /2ir 

=  (T^yifT  )/1.  108  (11) 

The  roll  damping  coefficient  is  derived  from  the 
observation  that  the  damping  ratio  of  a  typical  ship  is 
approximately 

r  =  0.1/ir  =  0.0318  (12) 

ship 

3y  working  backwards  through  the  second-order  system 
relationships,  the  roll  damping  coefficient  was  found  to  be 


2 

B  =0.1  (I.  A  G 23  )  /  ir  (13) 

S’ 

On  the  forcing  function  side  of  the  equation  of  motion, 
following  the  classical  usage  of  Froude  and  Krylov, 
[^Chadwick  and  Klotter,  1953J  the  ship-sea  coupling 
coefficient  was  considered  to  be  approximately  equal  to  the 
static  righting  coefficient.  Therefore,  the  equation  of 
motion  for  a  ship  in  a  seaway  is: 


J  i'  ♦  B  i  ♦  K  0  =  K  4>  (14) 

s  s  s  s 

Dividing  through  by  the  roll  inertia  coefficient  and 
performing  the  Laplace  Transform  operation,  the  system 
transfer  function  was  obtained: 


(B  /J  )  s  ♦  (K  /J  )1  <J»  ( s )  =  (K  /J  )  vp(s) 
s  S  s  S  -J  s  s 

2  _  2 

(2  f  w  )  s  ♦  «  1  <j>  (S)  =  w  q,  (S) 

S  S  S  J  s 


(15) 

(16) 


*  (s) 


*  (s) 


~2 
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ui 

_ s _ 

♦  |2t«)s 
s  s 
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~2 

iii 
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(K  /J  ) 


(17) 


s  ♦  (B  /J  )  s  «•  (K  /J  ) 
s  s  s  «. 

It  was  useful  in  the  analysis  to  normalize  the  equations 
such  that  the  ship's  natural  frequency  was  unity.  This  was 
accomplished  in  the  following  manner: 

Starting  with  the  general  equation  of  motion,  divide 
through  by  the  static  righting  coefficient: 

(J  /K  )  0  ♦  (B  /K  )  <p  +  q>  -  4*  (18) 

s  s  s  s 

Let  t'  =  a.  t  =  (/X  /J  )  t 
s  s  s 

Def  ine  d  0  /dt  =  <t> 

d  <t>  /dt '  =  0' 

2  2 
d  0  /dt '  =  <0 

Then  <t>  *  (JK  /J  )  <t>‘ 
s  s 

0  =  (K  /J  )  <*>» 
s  s 

The  ejuation  was  then  rewritten  as 

0"  ♦  (B  /  K  J  )  <0*  ♦  0  =  i  (19) 

s  s  s 

The  LaPlace  Transform  operation  was  performed,  and  the 
system  transfer  function  (normalized)  was  ootained: 


[s2  ♦  (B  /  K  J  )s  *  1~]  *  (s)  =  *  (s) 

s  s  s  -J 


♦  (S) 

*  (S) 


s  ♦  (B  /  K  J  )  s  ♦  1 
s  s  s 


Frequency  Response. 


In  order  to  study  the  frequency  response  of  a  vessel 
to  the  waveslope  spectrum  of  the  ocean,  it  was  necessary  to 
introduce  a  numerical  example.  The  example  chosen  was  a 
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well  used  example  ^Blagoveshchensky ,  1962J,  which  was: 

L  =  49.1  meters  (161.0  feet) 

b  =  9.0  meters  (  29.5  feet) 

A  =  936.0  metric  tons  (921.2  long  tons) 

GM  =  0.73  meter  (  2.4  feet) 

=  3.63  seconds 

The  radius  of  gyration  of  the  mass  of  the  ship,  and  the 
coefficients  of  the  equation  of  motion  were  computed  as 
follows : 

k  =  (8 . 63 v/( 9.  81)  (0.73)  )  /2*  =  3.63  meters 

2  2 
J  =  (936/9.  81 )  (3 . 63)  =  1289.021  tons-meters-sec 

s 

2 

B  =  (0.  1)  (3.63)  (936)  (0.73)/*  =59.7  tons-meters-sec 

3 

K  =  (936)  (0.7  3)  =  683.28  tons-meters 
s 

These  values  were  then  substituted  into  the  transfer 
function,  the  roots  of  the  characteristic  equation  were 
found,  and  the  frequency  response  was  determined. 


4*  (s)  =  _ ^ _ 0.53 _ 

q,  (s)  s  ♦  0.046s  ♦  0.53 


(22) 


roots  (poles)  at  s  =  -0.02  ±  j  0.73 

natural  frequency  ui  =  0.73 

s 

damping  ratio  f  =  0.032 

s 

The  frequency  response  is  graphically  illustrated  in 
Figure  1-2.  The  Digital  Simulation  Language  (DSL/360) 
program  which  generated  the  plot  is  listed  in  Appendix  A. 
The  zero  decibel  low  frequency  response  demonstrates  the 
natural  tendency  of  a  ship  to  ride  normal  to  the  surface  of 
a  wave.  The  resonant  peak  which  occurs  at  0.724  rad/sec  has 
a  gain  of  23.9  db,  or  a  magnification  factor  of  15.6. 
According  to  the  waveslope  spectrum  of  the  ocean,  tabulated 
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FREQUENCY  RESPONSE ,  SHIP  WITHC'J 

Closed  loop  response,  nqg lCB*  v 


Figure  1-2.  Frequency  Response  of  Typica  Ship 


in  Appendix  E,  for  a  sea  of  15  foot  significant  wave  height 
(sea  stats  6)  the  maximum  waveslope  angle  at  0.724  raa/sec 
would  be  2.23  degrees.  Therefore,  the  ship  would  experience 
a  maximum  roll  angle  of  34.3  degrees,  relative  ro  the 
horizon.  This  severity  of  roll  causes  crew  discomfort,  and 
limits  the  ability  of  the  ship  to  conduct  operations,  such 
as  gunfire  support,  flight  operations,  or  small  boat 
operations. 
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11  •  SHIP  STABILIZATION  BY  PASSIVE  ANTI-ROLL  TANKS 

A.  DEVELOPMENT  OP  THE  EQUATIONS  OF  MOTION 

In  order  to  apply  servo  theory  to  the  problem  of  ship 
stabilization#  the  physical  system  was.  assumed  to  be  linear, 
with  lumped  parameters  and  constant  coefficients.  It  was 
also  assumed  that  the  rolling  degree  of  freedom  could  be 
considered  independently  of  the  other  degrees  of  freedom, 
such  that  the  ship  oscillates  as  if  about  a  fixed  center  of 
rotation.  The  problem  thus  became  a  two-body  problem,  the 
ship  and  the  tank-fluid,  each  with  one  degree  of  freedom. 
The  tank-fluid  was  idealized  as  a  one-dimensional  filament, 
with  provision  for  varying  cross-section  but  with  constant 
velocity  across  each  cross-section.  These  assumptions 
required  small  angles  of  ship  roll  but  allowed  for  moderate 
angles  of  tank-fluid  level,  both  being  reasonable  in  a 
stabilized  ship. 

1 .  Doable  Pendulum  Analogy. 

As  discussed  earlier  in  this  paper,  a  ship  at  sea 
behaves  like  a  damped  oscillator,  and  was  modeled  as  a 
second-order  spring-mass-damper  system.  Another  appropriate 
model  which  was  considered  was  the  pendulum  model  ^Chadwick, 
1950j  [^Chadwick  and  Klotter,  1953^j,  with  a  spring  constant, 
an  inertia  constant,  and  a  damping  constant. 

The  stabilizer  chosen  was  a  U-tube  anti-rolling  tank 
system.  Tho  water  or  other  fluid  in  the  U-tube  acts  as  a 
hydraulic  pendulum.  Since  the  tank  was  to  be  contained 
within  and  wholly  supported  by  the  ship,  the  ship-tank 
system  that  resulted  was  a  pendulum  supporting  a  pendulum, 
or  a  double  pendulum.  In  general,  this  kind  of  system  leads 
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to  two  second-order  differential  equations  with  coupling 
ter  os . 


2 .  Noiaencla  tura. 

There  are  two  types  of  nomenclature  needed  here. 
These  are  the  parameters  and  variables  of  the  differential 
aquations,  and  the  dummy  variables  and  parameters  which 
characterize  the  tank  geometry. 

a.  Differential  Equations. 

The  differential  equations  are  second-order  in 
nature,  and  have  forcing  functions.  The  variables  and 
parameters  are  as  follows: 

Dependent  variables  of  the  forcing  side: 

4>  -  instantaneous  effective  waveslope 
Parameters  of  the  forcing  variables: 

K  =  ship-sea  coupling  coefficient 

S3 

Torques  produced  on  the  forcing  side: 

K  ^  =  disturbance  torque  due  to  waves 

ss 

Dependant  variaoles  of  the  homogeneous  side: 

<t>  -  ship  roll  angle  relative  to  the  horizon 
0  =  tank-fluid  level  angle  relative  to  the  ship 
Parameters  of  the  homogeneous  variables: 

Ship  self -terms: 

J>  =  roll  inertia  coef.,  unmodified  ship  (no  tanks) 

3 

B'  =  damping  coef.  of  unmodified  ship 

3 

K*  =  static  righting  coef.  of  unmodified  ship 

3 

J  =  roll  inertia  coef.,  modified  ship  (with  tanks) 

3 

3  =  damping  coef.  of  modified  ship 

3 

K  =  static  righting  coef.  of  modified  ship 
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Tank  self-tarns : 


J  =  equivalent  inertia  coaf.  of  tank-fluid 
=  skin-friction  damping  coef.  of  tank-fluid 
3  =  overall  damping  coaf.  of  tank-fluid 

K  =  equivalent  static  coef.  of  tank-fluid 

Mutual  coupling  terms: 

J  =  mutual-coupling  inertia  coefficient 
5  t 

3  =  mutual-coupling  damping  coefficient 

3 1 

K  =  mutual-coupling  static  coefficient 

3  t 

It  is  a  feature  of  the  double  pendulum  that 

X  =  K 
3 1  t 

B  =  0 
St 

As  discussed  earlier  in  this  paper 

K  =  K 
3  3  S 

b.  Tank  Geometry. 

The  notation  used  in  the  derivation  whicn 
depends  upon  tank  geometry  is  shown  in  Figure  2-1,  and  is 
given  below. 

s  «  a  generalized  coordinate,  the  curvilinear 

distance  along  the  fluid  trajectory.  s=0  at 
the  tank-fluid  surface  in  the  starboard  tank, 
and  s=S  at  the  tank-fluid  surface  in  the  port 

O 

tank,  when  9  =0 . 

H  =  lever  arm  of  tank,  the  horizontal  distance 

from  the  vertical  centerline  of  either  side 
tank  to  the  vertical  centerline  of  the  ship. 
r(s)  =  the  radial  distance  from  the  virtual  center 
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Figure  2.1.  The  Dependent  Variables  and  Nomenclature  of  Tank 
Geometry 
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of  rotation  to  any  point  on  the  fluid 
trajectory. 

P  =  density  of  fluid  in  tanks. 

D(s)  =  the  perpendicular  distance  from  the  virtual 
center  of  rotation  to  tangent  at  any  point  on 
trajectory. 

a(s)  =  angle  between  r(s)  and  D  (s)  . 

•y  ( s)  =  angle  between  r  (s)  and  deck  of  ship. 

A  =  tank  cross-section  at  s=0  and  s=S. 

a(s)  =  tank  cross-section  at  any  point  s. 

g  =  acceleration  of  gravity. 

3.  Euler-LaGrange  Equations  of  Motion. 

Tie  two  degree  of  freedom  system  of  ship  and 
tank-fluid  was  best  treated  by  means  of  the  Euler-LaGrange 
equations  in  generalized  coordinates.  The  generalized 
aquation  is  [jchadwick  and  Klotter,  1953^ 

d  rail  a  t  a  p  a  p 

dt  aq  ♦  aq  ♦  aq  +  a  q  =  Q  (23) 

ij  i  i  i  i 

where  is  the  so-called  generalized  ncn-conservati ve  force 
in  the  ith  equation,  and  i«1,2.  The  generalized  force  is 
equal  to  the  rate  of  doing  virtual  work  against  the  external 
non-conservative  forces  when  q^  is  given  a  virtual 
displacement.  The  parameters  of  the  Euler-LaGrange  equation 
are : 


T  = 

2 

Z  mv  /2 

=  kinetic  energy 

(24) 

0  = 

z  ftmix 

=  conservative  potential 
(due  to  gravity) 

energy 

(25) 

D  = 

Z  0  v2/2 

=  dissipative  energy 
(due  to  skin  friction. 

,  4  (26) 
etc. ) 

f  = 

conservative 

force 

0  =  damping  parameter 
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Calculation  of  the  Energy  Functions. 


For  each  -ype  of  energy  (kinetic,  conservative 
potential,  and  dissipative)  ,  the  contribution  of  the  ship 
alone  and  the  contribution  of  the  tank-fluid  alone  was 
calculated,  then  the  contributions  were  added  to  arrive  at 
the  total  energy.  In  the  regions  where  the  fluid  surfaces 
travel,  namely  the  vertical  legs  of  the  (J-tube,  it  was 
assumed  that  the  tanks  were  straight,  parallel,  and  of 
uniform  cross-section.  This  meant  that  a (s)  and  D(s)  were 
constants  in  these  regions,  thus  for  any  expressions  which 
depended  only  on  the  variables  a(s)  and  D(s),  the  limits  of 
integration  were  fixed  at  zero  and  S,  without  approximation. 

Calculation  of  kinetic  energy: 


fluid 


S+H  tan d 


u  *  /"-4^{p*TtT]2,[ir(s)]2*  2p^rtr][iD<s>] 


R  tans 


ship 


total 


Jf  <*>  /2 
s 

r  +  t 

fluid  ship 


Calculation  of  potential  energy: 


5 

°fluid  ~  p  3  A  R  y  a^s)  r  js)  sin7  ds  1  -  cos  <t> 

L  A  J  L  - 


2  r  2 

♦  p  g  A  R  2  tan  0  sin  <t>-  tan  d  cos  <t> 


2  L 
U  =  K'  $  /2 

ship  s 

a  =  o  ♦  a 

total  fluid  ship 
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Calculation  of  dissipative  energy: 


D 

fluid 

.2 

B«  0  /2 
t 

(33) 

D  = 

•2 

B’<*>  /2 

(34) 

ship 

s 

D 

D  ♦  D 

(35) 

total 

fluid  ship 

5.  The  Equations  of  .Motion. 

The  equations  of  motion  of  the  ship  and  tank,  system 
come  from  equation  (23)  ,  with  qj  =  <(>  ,  and  =  9  . 

Substituting  in  equations  (27)  through  (35)  and  simplifying, 
and  letting  Q  =  K  t  and  Q,  =  0,  the  equations  of  motion  are 

Shin:  (J  *  *3  <j>  *K  <t>  )  +  (J  0  *K  9  )  -  K  *  (36) 

s  s  s  St  t  s 

Tank:  (J  *  *K  <t>  )  *  (J  9  *B  9  +K  9  )  *  0  (37) 

st  t  t  t  t 

Obtaining  the  system  transfer  function  from  the 
equations  of  motion  was  somewhat  more  complicated  than  it 
was  for  the  ship  without  tanks.  First,  the  equations  were 
divided  through  by  the  respective  inertia  coeficients,  and 
the  LaPlace  Transform  operation  was  performed: 

2  2 

Ship:  fs  ♦  (B  /J  )  s*  (K  /J  fU(s>  +  f(J  /J  )s  ♦  (K  /J  fle(s) 
*-  ss  ss-1  '-sts  t  s  -* 

=  (K  /J  )  *(S)  (38) 

s  s 

2  2 

Tank:  [~(J  /J  )s  ♦  (K  /J  )]<t>  (s)  t  fs  ♦  (B  /J  )  s*  (K  /J  )"le(s) 
Lstt  ttJ  L  tt  tt-1 

■  0  (39) 

$  (s)  was  isolated  by  applying  Cramer’s  Rule,  and  the 
transfer  function  was  obtained: 
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4 

B  J  +J  B 

St  St 

3 

r  1 

K  J  +J  K  +B  B  -2K  J 

=  s  «• 

l 

s  ♦ 

St  St  jj  t  tst 

J  J  -J 

J  J  -  J 

St  St  _ 

_  S  t  St 

B  K  +K  B 
St  St 

r  2  ~1 
K  K  -K 
s  t  t 

2 

J  J  -J 

St  St 

s  ♦ 

2 

J  J  -0 
_  s  t  st_ 

2 

/ 3  )  s  +  (B 

/J 

)  s+ (K  /J  ) 

s  s 

t  z 

t  t 

*  (3) 


2 

s 


(40) 


(41) 


6.  Definitions  of  Parameters. 


The  parameters  of  the  equations  of  motion  are  given 
in  terms  of  previously  define!  parameters  and  dummy 
variables  as  follows: 


Ship  self-terms: 

S 

Js  =  j'  ♦  P AH  y"[ a  (s)/a]  Q:  (s) /R^  ds 
0 


B  =  B* 
s  s 

S 

Kg  a  K*  ♦  pgAH  /[a(s)/A][r(s,/H]  sin  y 

o 

Tank  self-terms: 

s 

J  =  P  A  R  J  A/a  (s)  ds 
0 

B  =  B' 
t  t 

2 

K  =  2  P  g  A  R 

t 


ds 


(42) 

(43) 

(44) 

(45) 

(46) 

(47) 
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Mutual-coupling  terms: 

S 

2  f 

(48) 

J  =  p  A  R  /  D  (s)  /R  ds 

st  J 

0 

B  =  0 

3 1 

(49) 

K  =  K 

3 1  t 

(50) 

B.  ANALYSIS  OF  THE  EQUATIONS  OF  MOTION. 

Having  arrived  at  a  set  of  equations  of  motion,  it  was 
necessary  to  conduct  a  brief  analysis  of  the  equations  to 
give  them  physical  significance. 

1 .  The  Self  Terms. 

Looking  first  at  the  self-terms  of  the  homogeneous 
equations,  it  was  observed  that  the  ship  acts  as  a  simple 
pendulum  if  the  tank-fluid  does  not  move,  and  that  the 
tank-fluid  acts  as  a  simple  pendulum  if  the  ship  is 
prevented  from  moving.  Each  of  these  pendulums  has  its  own 
natural  frequency.  These  two  frequencies  are  important 
factors  in  the  stabilization  problem. 

2.  The  Mutual  Terms. 

There  is  a  third  critical  frequency  associated  with 
the  ship-tank  system,  called  the  secondary  resonance 
frequency,  which  is  due  to  the  mutual-coupling  terms.  The 
mathematical  significance  of  the  effect  became  apparent  when 
the  ship  aquation,  equation  (36) ,  was  rewritten  with  the 
mutual  teems  on  the  right  side  as  follows: 

J  i>'*B  6  J  *K  J)  (51) 

3  s  s  s  St  t 
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This  aquation  shows  the  ship  as  a  simple  pendulum,  acted 
upon  by  a  net  torque,  which  is  the  torque  due  to  waves  minus 
the  torque  due  to  tank-fluid  motion.  If  the  tank-fluid  were 
made  to  oscillate  at  the  secondary  resonance  frequency,  the 
tank-fluid  torque  would  be  equal  to  the  wave  torque,  and 
thus  the  ship  would  be  influenced  by  zero  net  torque. 


C.  ANT I- ROLL  TANK  DESIGN  CRITERIA. 


The  easiest  tank  geometry  to  design  and  build  is  a  right 
rectangular  0-tube,  with  the  vertical  legs  having  constant 
cross-sectional  area,  and  the  horizontal  cross-leg  having 
constant  but  smaller  cross-sectional  area.  Given  this  tank 
geometry. 


a(vertical  leg) 
a  (horizontal  leg) 

D(vertical  leg) 

□  (horizontal  leg) 


A 

a 

2  2 

(A  R  w  )/(g  -  h  w  ) 
s  s 

R 

D 


(52) 

(53) 

(54) 

(55) 


The  length  of  the  trajectory  of  the  fluid  in  the  tank, 
S,  is  dependent  upon  the  lever  arm  of  the  tank,  R,  and  the 
height  of  the  fluid  column  in  a  vertical  leg  of  the  O-tube 
when  9=0°,  h, 

S  =  2  (R  ♦  h)  (56) 

For  maximum  effect,  the  tank  system  must  span  the 
extreme  beam  of  the  ship.  Therefore  the  lever  arm  of  the 
tank  is  half  of  the  beam,  b,  diminished  by  the  width  of  a 
vertical  leg,  w, 

R  =  (b-w)  /2  (57) 

The  relationships  derived  from  the  equations  of  motion 
which  define  the  tank  dimensions  and  parameters  are: 
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K  /J 
t  s 


X 


(58) 


VJt 


2 

CJ 

t 


/ 


_ l  .3- 

A/a  <s) 


ds 


J  /J  =  X 
st  t 


X 


=  J  (D  a)  /  (R  A)  ds 
0 

=  2  h  *  2  ( D  a/A) 


(59) 


(60) 


The  first  ratio.  A,  reflects  the  free  surface  effect  of 
the  tank-fluid,  and  should  be  a  small  number.  An  acceptable 
value  was  found  to  be  0.05  [jfachanaratana ,  1973J. 

To  be  effective  as  a  stabilizer,  the  anti-roll  tank 
system  must  be  "tuned”  to  the  natural  frequency  of  the  ship 
in  which  it  is  installed.  Therefore, 


2  2 

U>  —  U1 

t  S 

K  /J  =  K  /J  (61) 

t  t  s  s 

The  third  ratio,  X,  is  dependent  upon  the  placement  of 
the  tanks  in  the  ship  relative  to  the  center  of  rotation. 
This  ratio  can  be  forced  to  zero  (J  forced  to  zero)  by 
placing  the  horizontal  leg  of  the  U-tube  above  the  center  of 
rotation  by  the  amount 


-D  =  h  A/a 


(62) 


There  are  no  hard  and  fast  rules  for  determining  a  value 


33 


of  h,  except  that  the  tanks  must  fit  within  the  hull  of  the 
ship,  and  the  fluid  in  the  tanks  should  not  be  permitted  to 
either  contact  the  top  or  completely  drain  out  of  the  bottom 
of  a  vertical  leg  of  the  tank  at  the  maximum  expected 
tank-fluid  level  angle,  ^  a  j  . 

Since  J  depends  upon  tank  geometry,  the  approximation 
must  be  made  that  =  J'  ,  for  purposes  of  initiating 

design.  Then 

K  =  X  J» 
t  s 

2 

=  2  p  g  A  R 


2 

=  P  g  1  w  (b-v)  /2  (63) 

A  reasonable  value  of  tank  width  should  be  chosen.  Then 


2 

1  =  2K^_/pgw  (b-v)  (64) 

flaking  the  approximation  that  the  natural  frequency  of 
the  modified  ship  is  the  same  as  the  natural  frequency  of 
the  ship  without  tanks,  and  picking  a  value  for  h,  equation 
(53)  is  solved  for  a,  and  then  aquation  (62)  is  tested  for 
feasibility.  If  the  value  of  D  so  computed  cannot  be 
physically  realized,  then  a  realizable  value  of  D  must  be 
chosen.  This  merely  means  that  JS|  will  not  be  identically 
zero,  and  the  overall  analysis  will  take  more  care. 

The  tank  damping  parameter  is  derived  by  working 
backwards  through  the  second-order  relationships  for  the 
hydraulic  pendulum,  noting  that  the  damping  ratio  for  the 
tank-fluid  is  approximately  equal  to  the  damping  ratio  of 
the  (unmodified)  ship: 


f 


t 


0 .  1/  w 


(65) 


B  =  0.2  ,/K  J  /  i r  (66) 

t  v  t  t 

Once  the  tank  system  parameters  have  been  determined, 
the  modified  ship  parameters  are  found.  Over  the  range  s=0 
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I 


I 


over  the 


The  design  steps  should  then  be  retraced,  removing  the 
approximations.  Normally,  after  two  or  three  iterations  of 
the  design  procedure,  the  calculated  values  converge  on  the 
actual  values,  and  a  final  design  has  been  achieved.  The 
FORTRAN  subroutine  used  in  this  iterative  design  procedure 
is  included  in  Appendix  B  to  this  thesis. 


D.  FREQUENCY  RESPONSE  OF  PASSIVE  SHIP-TANK  SYSTEM. 

The  same  numerical  example  used  in  the  frequency 
response  study  of  a  vessel  on  the  ocean  was  used  for  the 
design  example  and  frequency  response  study  of  a  ship  with 
passive  anti-roll  tanks  to  the  waveslope  spectrum  of  the 
ocean : 

2 

J*  3  1239.021  metric  tons-met ers-sec 

3 

B*  3  59.746  metric  tons- meters-sec 

3 

K'  3  683.280  metric  tons-meters 

3 

The  following  arbitrary  choices  were  made: 
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=  i.O  i  e  t  <3  r 
'  .  3  a  •9*-  er 
=  i .  Jc 

3 

d  :  -  i  .  jlo  .litri;  tons/meter  (salt  water  used) 

1  .  «-*  ~ 

'f  =  (0.05)  (1  283.021) 

=  64.431  metric  tons-aeter 

2 

1  =  (2)  164.  45  1)  /  (1 .025)  (1)  (9-  1) 

=  1.965  aeter 

A  =  w  1 

2 

=  1.965  aeter 

2 

=  0.457  aeter 

Testing  equation  (52)  for  feasibility, 

-D  =  (1.3)  (1.965)  /  0.457 

=  5.59  aeter 

Since  this  tank  placement  would  probably  not  be  possible 
in  the  ship  being  studied,  a  more  realistic  value,  still 
above  the  center  of  rotation,  was  arbitrarily  chosen: 

-0  =  1.0  Deter 

Then  X  is  non-zero,  as  indicated. 

J  /J  =  X 
3 1  t 

=  2  [j.  3  -  (0.457)  (1.965)] 

=  2.135 

To  suaiarize,  the  tank  dimensions  found  thus  far  were 
w  =  1.0  meter 

1  =  1.965  meter 
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A 


2 

=  1.965  meter 

2 

a  =  0.457  meter 

H  =  4.0  meter 

D  =  -1.0  meter 

Siven  these  tank  dimensions,  the  modified  ship 
parameters  became: 

J  =  1289.021  ♦  11.155 

5 

2 

=  1300.176  metric  toQs-meters-sec 

K  =  683.280  -  3.121 

3 

=  680.159  metric  tons-meter 


After  four  iterations  of  the  design  procedure,  the 
following  final  values  were  obtained: 
w  =  1.0  meter 

1  =  1.982  meter 

2 

A  =  1.982  meter 

2 

a  =  0.454  meter 

R  =  4.0  meter 

D  =  1.0  meter 

2 

J  =  1300.236  metric  tons-met ers-sec 

3 

B  =  59.746  metric  tons-met ers-sec 


K  ~  680.135  metric  tons-raeters 

3 

2 

J  =  124.276  metric  tons-meters-sec 

t 

B  *  5.722  metric  tons-meters-sec 

t 

K  =  65.012  metric  tons-meters 

t 

2 

J  =  266.152  metric  tons-meters-sec 

3  t 

These  values  ware  then  substituted  into  the  transfer 
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function,  the  roots  were  located,  and  the  frequency  response 
was  obtained. 


2 

*  (s) 

*  (s)  s  *  0.164  s  +  1.485  s  ♦  0.0857  s  +  0.441 

(73) 

zeroes  at  s  =  -0.0230  t  j  0.723 

pales  at  s  =  -0.0135  t  j  0.642 

and  s  =  -0.0684  t  j  1.031 

The  frequency  response  is  graphically  illustrated  in 
Figure  2-2.  The  DSL/360  program  which  generated  the  plot  is 
listed  in  Appendix  3.  The  negative  decibel  low  frequency 
response  demonstrates  the  ability  of  the  passive  tank  system 
to  overcome  the  natural  tendency  of  the  ship  to  ride  normal 
to  the  waveslope,  making  it  ride  upright  with  respect  to 
earth  reference  instead. 

Comparing  this  response  with  the  response  of  the 
unmodified  ship,  given  in  Figure  1-2,  it  was  observed  that 
the  effect  of  adding  passive  anti-roll  tanks  was  to  place  a 
pair  of  complex  zeroes  at  the  ship's  own  natural  frequency 
to  reduce  the  resonant  peak,  from  24  db  at  0.7  rad/sec,  to 
two  lower  peaks  of  14  db  at  0.6  rad/sec  and  10  db  at  1.0 
rad/sec,  and  a  notch  of  - 1 1  db  at  0.7  rad/sec.  In  sea  state 
6,  this  translates  to  maximum  ship  roll  angles  of  8.5 
degrees  at  0.6  rad/sec,  and  15.3  degrees  at  1.0  rad/sec,  a 
considerable  improvement  over  the  34.8  degree  maximum  roll 
angle  of  the  unmodified  ship. 
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FREQUENCY  RESPONSE,  PASSIVE  TfiNKS 
§Cl_OSED  LOOP  RESPONSE,  MPG  !DB)  VS.  LOG 


III.  SHIP  STABILIZATION  SI 


ACTIVATED  ANTI-ROLL  TANKS 


A.  MODIFIED  EQUATIONS  OF  MOTION. 

Activiting  an  anti-roll  tank  system,  n  ithemar ically , 
only  involves  aiding  a  forcing  function  ters  to  the  tank 
equation  of  motion.  Physically,  hrwever,  it  can  be 
accomplished  in  basically  two  different  ways.  These  methods 
are  pneumatic  and  hydraulic.  A  pneumatic  effector  operates 
by  alternately  creating  an  overpressure  or  a  partial  vacuum 
condition  in  the  airspace  above  the  tank-fluid  surface  in 
the  vertical  legs  of  the  U-tube,  while  an  hydraulic  effector 
generally  consists  of  an  arrangement  of  pumps  in  the 
horizontal  cross  duct.  The  effector  type  chosen  for  this 
analysis  was  a  constant  speed  controllable  pitch  propellor 
p-*-2,  so  chosen  for  its  linear  blade  pitch  angle  vs.  flow 
characteristic  [^Blagoveshchensky,  1962^. 

1 .  Nomenclature. 

In  addition  to  the  nomenclature  of  the  passive  tank 
system,  the  following  are  the  parameters  and  variables  of  an 
activated  tank  system: 

Dependent  variables  of  the  forcing  side: 

a  -  instantaneous  pump-blade  pitch  angle  (commanded) 

Parameters  of  the  forcing  variable: 

B  =  velocity  head-loss  coefficient  of  the  pump 
P 

K  =  static  head  coefficient  of  the  pump 

3 

Torques  produced  on  the  forcing  side: 

K  a  -  B  9  -  Counter-rolling  torque  due  to  pump 

?  P 
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2 .  Modification  of  Parameters. 

As  indicated  by  the  expression  for  counter-rolling 
torque,  the  addition  of  a  pump  to  the  tank  system  has  the 
effect  of  increasing  the  overall  damping  coefficient  of  the 
tank-fluid.  All  other  parameters  of  the  homogeneous  side 

remain  unchanged.  The  pump  should  be  chosen  such  that  the 
composite  tank  damping  factor  is  0.45  [j/achanaratana,  1973J: 

B  =  B*  *  B  (74) 

t  t  p 

B  /J  »  0.45  (75) 

t  t 

3  =  0.45  J  -  3*  (76) 

p  t  t 

3.  Equations  of  Motion. 

The  equations  of  motion  for  the  activated  anti-roll 

>2  *  K 

(77) 

S 

Tank:  (J  0  +  K  <t>  )  *  (J  M  B  J  +  K  9  )  =  Ka  (78) 

st  t  t  t  t  p 

Following  the  procedure  used  for  the  passive  ship-tank 


(23) 

,  with 

*  Ki* 

and 

Klotter , 

1953]' 

(J 

y  ♦  k 

8  ) 

-  K  4> 

st 

t 

s 

♦  B 

9  +  K 

8  ) 

=  K  a 

4»  (s) 

was  found 

to  be 

<t>  (s) 

=  *  (s)  - 

*  (s) 

(79) 

+  (s) 

=  (K  /J  ) 
s  s 

fs2  ♦  (B  /J  )  s  ♦  (K  /J  )1  *  (s) 

L  t  t  t  t  -J 

(80) 

A 

*2(s> 

=  (K  /J  ) 
P _ - 

f(J  /J  )  s2  ♦  (K  /J  )"]  A  (s) 

<-  st  s  t  s  J 

(81) 

A 

a 

=  as  defined  by  equation  (40) 

B.  CONTROLLER  DESIGN  AND  FREQUENCY  RESPONSE. 


The  problem  of  ship  roll  control  is  a  frequency  domain 
problem:  controlling  against  a  (time  domain)  stochastic 


disturbance  with  a  generally  known  (frequency  domain) 
spectrum.  Therefore  the  classical  approach  to  control  was 
pursued,  designing  on  the  basis  of  feedback  compensation. 
It  was  assumed  that  roll  angle,  rate,  and  acceleration 
signals  would  be  available  from  shipboard  sensors  (gyros  and 
accelerometers) .  The  controller  transfer  function  was 
considered  to  be  of  the  form 


C(s)  =  A(s)/+  (s) 


=  K  S  +  K  s  ♦  K 
2  1  0 


(82) 


For  purposes  of  block  diagram  development  and  analysis, 
the  following  block  transfer  functions  were  defined,  based 
upon  equations  (79),  (80),  and  (81|: 

SS(s)  =  4,^  (s)  /  qr(s) 

2 


=  K  /J  fs  ♦  (B  /J  )S  ♦  (K  /J  )~|/A 
s  s  ^  tt  tt~* 

2  r*  2  2-] 

=  w  s  ♦  2  f  «  s  +  «  /  i 

n  1-  t  n  nJ 


ST  (s)  =  d>  (s)  /  k  (s) 

2 


=  K  /J  r C J  /J  )s  ♦  (K  /J  )!/ 
pt'-sts  tsJ 


=  K  /J  [: 
p  s  *- 


2  21 
Xs  ♦  w  /  A 

n-1 


=  k  /k  r 

p  s 


2  2- 
Xs  *  w 


s  +  2  ( w  s  +  ui  1  (85) 

t  n  n  J 


SST  (s)  =  ST  (s)  /  SS  (s) 

:?c- 

K  /J  [  s*  (B  /J  )  s  ♦  (K  /J  )1 

s  s  ^ _  t  t _ t  t  J  - 

A  -  K  /J  fs  ♦  (B  /J  )  s  ♦  (K  /J  )1 

SS'-  tt  tt-* 


G(3) 


(83) 


(84) 


w  fs 
n  <- 

i 

♦ 

2  {*  cu  s  ^ 
t  n 

l — i 
a 
3 

2 

r  2 

2-1 

A  -  oj 

Is 

♦  2  f  <*»  s 

+  u> 

n 

u. 

t  n 

n  J 

(86) 
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H(s) 


Where 


=  SST(s)  C(s) 

r  4  3  222  ^  n 

=  1/K  XC  S  +  XC  S  <•  (XC  C  )  S  ♦  w  C  S  +ui  C  1 
s1-  2  1  0  n  2  n  1  n  0-1 

2  2 
s  ♦  2  s  ^  s  ♦  w 
t  a  n 


-  K  K 
P  2 


2 

JL- 
(87) 


C  =  K  K 

1  p  1 

C  =  K  K 

3  p  0 

The  "open  loop"  transfer  function,  then,  was 
F(s)  =  G  (s)  /  [l  *  G  (s)  H  (s)] 


(K  /J  )  Ts  *  (B  /J  )  s  ♦  (K  /J  )~] 
s  s  L _ t  t _ t  t  -I 

A  (open) 


(88) 


A  ( open)  =  A  -  (K  /J  )  fs2  ♦  (B  /J  )  s  ♦  (K  /J  )~I  ♦ 

ss1-  tt  tt-1 

r  4  3 

♦  (1/J  )|XC  S  *X  C  s  ♦  (XC  ♦  (K  /J  )C  )s 
s  L  2  1  0  t  t  2 

>(K  /  J  )C  s  ♦(K/J  )C  “1 

t  t  1  t  t  oJ 


(89) 


The  overall,  or  "closed  loop"  response  became 

T(s)  =  <*>  (s)  /  *  (s) 

=  F(s)  /  [l  ♦  F(s)] 

*  (K  /J  )  fs2  ♦  (B  /J  )s  *  (K  /J  )1 

-  s  s  L _ t  t _ tt  J 

A  (closed) 


(90) 


r  4  3  2 

A  (closed)  =  A  «•  (1/J  )  jjC^s  ♦XC^  ♦  (XCq+ (K^ /J^)  C^)  s 

(91) 


♦  (K  /J  )  C  s  ♦  (K  /J  )  C  1 
t  t  1  t  t  0-1 
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The  block  diagram  development  is  illustrated  in  Pigure 
3-1.  The  first  design  approach  attempted  was  the  modern 
control  theory  state  variable  feedback  method.  The 
difficulty  with  this  method  was  the  inability  to  define  a 
time-domain  objective  or  cost  functional. 

The  next  design  approach  attempted  was  with  the  Computer 
Aided  Design  of  Optimal  Compensators  (CADOC)  program 
[  Anastassakis ,  1977].  This  method  consisted  of  selecting  a 
frequency-domain  cost  functional  and  a  desired  frequency 
response  curve,  and  used  the  function  minimization  procedure 
of  M.J.  Box.  There  were  two  difficulties  encountered  in 
applying  this  program  to  the  problem  of  ship  anti-roll 
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(a)  Direct  implementation  of  equations  (79)  and  (82). 


(b)  First  state  of  block  diagram  reduction. 


Figure  3.1.  Block  Diagram  Development  for  Activated 
Ship/Tank  System 


(d)  Equivalent  block  diagram  for  servo  theory  analysis 


Figure  3.1. 


Block  Diagram  Development  for  Activated 
Ship/Tank  System 


control.  The  first  problem  was  that  a  portion  of  the 
feedback  path  in  the  physical  system  was  fixed  by 
architecture,  while  the  CADOC  program  required  all 
coefficients  in  the  feedback  path  to  be  adjustable.  The 
second  problem  involved  selecting  an  attainable  frequency 
response  for  CADOC  input.  It  was  concluded  that  the  CADOC 
program  was  best  suited  for  optimizing  a  close  design, 
rather  than  for  initiating  a  design. 

Another  design  approach  attempted  was  in  the  root  locus 
parameter  space.  Two  separate  procedures  exist  here,  and 
both  were  attempted.  The  first  procedure  was  the  standard 
root  locus  approach,  fixing  two  of  the  three  adjustable 
gains  and  varying  the  third.  This  method  amounted  to  an 
exhaustive  search  technique,  attempting  to  locate  the  best 
combination  of  the  three  gains  to  define  system  roots  with 
adequate  damping. 

The  second  root  locus  parameter  space  procedure  involved 
using  a  parameter  plane  computer  program,  PARAM  A,  which 
exists  at  the  Naval  Postgraduate  School  and  is  available 
from  Dr.  3.J.  Thaler,  but  has  not  been  documented  in  any 
literature  to  his  knowledge.  This  program  varies  two  gains 
which  appear  in  linear  combinations  in  the  coefficients  of 
tha  system  characteristic  equation,  for  a  fixed  value  of  the 
third  gain.  The  output  of  this  program  is  a  two-dimensional 
parameter  plane  showing  families  of  curves  of  constant 
damping  ratio  and  constant  natural  frequency  as  functions  of 
the  two  gains. 

The  problems  associated  with  both  root  locus  parameter 
space  procedures  were  the  exhaustive  search  nature  of  a 
three-variable  system,  and  the  exclusions  of  the  system 
closed-loop  zeroes  from  the  analysis.  This  was  noted  to  be 
a  problem  when  the  gains  found  by  these  procedures  were  used 
in  a  frequency  response  study  of  the  activated  ship/tank 
system,  and  large  resonant  peaks  were  observed. 
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After  attempting  all  of  the  above  methods  and  abandoning 
them  for  the  reasons  stated,  an  exhaustive  search  of 
closed-loop  frequency  response  curves  was  conducted,  using 
the  program  listed  in  Appendix  C.  An  acceptable  design  was 
accomplisaed  when  C  =  110,  C  =  -150,  and  C  =  -80.  This 
corresponded  to  placing  the  compensator  roots  at  s  =  -2.4388 
and  s  -  *0.5638,  which  improved  the  overall  system  damping 
and  provided  the  frequency  response  illustrated  in  Figure 
3-2.  This  response  was  characterized  by  two  low  resonant 
peaks,  one  at  0.582  rad/sec  of  2.85  db  (linear  gain  1.39), 
and  one  at  1.38  rad/sec  of  1.96  db  (linear  gain  1.25).  The 
minimum  response  between  the  two  resonances  occured  at  0.861 
rad/sec,  and  was  -4.75  db  (linear  gain  0.58). 

In  sea  state  6,  the  maximum  waveslopes  at  the  resonant 
frequencies  are  1.44  degrees  at  0.582  rad/sec,  and  8.11 
degrees  at  1.38  rad/sec,  and  the  maximum  waveslope  at  0.861 
rad/sec  is  3.15  degrees.  Using  this  controller  design,  the 
maximum  roll  angles  that  the  ship  should  experience  in  sea 
state  6  would  be  2.0  degrees  at  0.582  rad/sec,  1.83  degrees 
at  0.861  rad/sec,  and  10.1  degrees  at  1.38  rad/sec.  The 
wave  height  power  spectral  density  at  those  frequencies, 
according  to  the  spectra  of  Pierson  and  Moskowitz,  are  205 
ft2  /Hz  at  0.582  rad/sec  (0.093  Hz),  80  ft  2  /Hz  at  0.861 
rad/sec  (0.137  Hz),  and  12  ft2/Hz  at  1.38  rad/sec  (0.219 
Hz) .  Therefore  it  was  concluded  that  even  though  the  pure 
numbers  indicated  a  poor  roll  control  at  the  higher 
frequency  resonant  peak,  the  waves  at  that  frequency  lack 
sufficient  power  to  cause  the  ship  to  attain  the  maximum 
theoretical  roll  angle. 
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Figure  3-2.  Frequency  Response  of  Activated  Ship/Tank 

System,  Proportional  P^us  Rate  Plus  Acceleration 
Control 
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IV.  EFFECTS  0?  VESSEL  DISPLACEMENT  VARIATIONS 

The  analysis  of  ship  motions  and  roll  control  thus  far 
had  been  conducted  under  the  assumption  that  the 
displacement  of  the  vessel  was  constant.  Over  the  short 
term,  the  assumption  may  be  true  enough,  but  for  a  ship 
underway,  consuming  fuel  and  stores,  the  assumption  cannot 
hold  true  beyond  a  few  days  at  best. 

As  a  ship's  displacement  decreases,  the  underwater 
surface  area  of  the  hull  decreases  also,  in  keeping  with  the 
Principle  of  Archimedes.  This  in  turn  leads  to  a  lowering 
of  the  Center  of  Buoyancy,  which  results  in  a  lowering  of 
the  Metacenter.  The  effect  on  the  metacentric  height, 
however,  is  not  so  clear,  as  the  exact  location  of  the 
Center  of  Gravity  in  the  ship  depends  more  on  the  placement 
of  the  onboard  weight  than  on  the  total  displacement.  In 
general,  fuel  tanks  are  located  low  in  a  ship;  thus 
consumption  of  fuel  is  removal  of  low  weight,  which  shifts 
the  Center  of  Gravity  upward,  decreasing  the  metacentric 
height.  It  should  be  noted  here  that  ballasting  empty  fuel 
tanks  with  sea  water  has  the  effect  of  restoring  the  ship  as 
much  as  passible  to  the  original  condition  of  statical 
stability  which  existed  at  the  beginning  of  the  voyage,  and 
extends  the  time  frame  of  applicability  of  the  assumption  of 
constant  displacement. 

Varying  the  ship's  displacement  therefore  varies  the 
natural  frequency  of  the  ship,  and  effectively  de-tunes  the 
anti-roll  tank  system.  The  only  physical  adjustment  that 
can  be  made  to  a  passive  tank  installation  underway  is  to 
vary  the  quantity  of  water  in  the  U-tube.  This  has  the 
effect  of  changing  the  limits  of  integration  of  the  design 
equations,  and  therefore  changes  the  parameters  of  the  tank. 
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Additional  adjustments  are  available  on  the  activated 
tank  system.  These  are  the  gain  of  the  pump,  which  is  a 
function  of  the  pump  rotational  speed,  and  the  feedback 
gains  from  the  sensors  to  the  pump  controls.  Thus,  if  a 
method  exists  for  determining  on-line  the  various  ship 
parameters  as  they  change,  the  controller  for  the  activated 
tank  system  can  be  made  adaptive.  Various  schemes  may  be 
considered,  from  Kalman  filter  type  parameter  estimators  to 
FFT  devices,  and  this  field  is  left  as  a  topic  for  further 


V.  CONCLUSIONS 


In  studying  ship  roil  control,  it  was  necessary  to 
become  familiar  with  the  nature  of  ocean  waves.  The 
waveslope  spectrum  of  the  ocean  which  was  generated,  when 
used  in  conjunction  with  the  Pierson-Moskowitz  wave  power 
spectra,  gave  much  insight  into  the  problem.  It  was  learned 
that  the  sea  surface  does  not  and  can  not  assume  the  slopes 
that  an  observer  on  a  ship  often  estimates. 

A  typical  ship  on  the  ocean  magnifies  the  waveslopes 
near  the  resonant  frequency  by  a  factor  of  about  15.  By 
installing  passive  anti-roll  tanks,  this  waveslope 
magnification  factor  was  split  into  two  resonant  peaks,  one 
with  a  factor  of  5,  and  the  other  with  a  factor  of  3. 
Activating  the  tank  system  further  decreased  the  waveslope 
magnification  factors  to  1.4  and  1.3,  respectively. 

Several  assumptions  have  been  made  over  the  years  in  the 
study  of  roll  control  by  (J-tube  tank  systems,  and  some  of 
them  were  sources  of  problems  in  the  design  of  physical 
systems.  The  first  assumption  was  that  a  ship  has  constant 
parameters:  displacement,  metacentric  height,  and  natural 
period.  This  assumption  holds  true  only  if  the  ship  does  not 
consume  any  fuel  or  stores,  or  if  these  weights  are  replaced 
with  appropriate  quantities  of  suitably  placed  ballast. 

Another  assumption  was  that  the  tanks  could  be  placed  in 
the  ship  in  such  fashion  as  to  drive  the  mutual -coupling 
inertia  coefficient  to  zero.  This  assumption  had  two 
limitations.  The  first  pitfall  was  the  direct  connection 
between  this  assumption  and  the  assumption  of  constant  ship 
parameters.  The  tecond  pitfall  was  that  the  tanks  had  to  be 
placed  in  an  impossiale  location  on  the  ship  to  accomplish 
the  objective  of  the  assumption. 
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Soil  control  by  pump  activated  tanks  is  perhaps  the  best 
form  of  roll  control  for  slow  moving  vessels,  or  for  vessels 
conducting  operations  either  at  slow  speed  or  stopped,  where 
activated  fin  stabilizers  lose  their  effectiveness.  It 
should  be  noted  that  a  pump  activated  tank  system  can  be 
used  to  cause  rolling,  as  an  aid  for  breaking  out  of  ice. 
It  can  also  be  used  to  correct  list,  including  damage-caused 
list,  as  a  safer  alternative  to  counter- flooding.  Therefore 
the  pump  activated  LJ-tube  tank  system  is  the  most  versatile 
form  of  coll  control  that  a  multi-purpose  vessel  can  have 
installed. 
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APPENDIX  A 


Program  Listing  for  Frequency  Response  of 
Unmodified  Ship,  with 
FORTRAN  Subroutine  PABAM 

The  Frequency  Response  program  was  written  in  Digital 
Simulation  Language  (DSL/360)  to  be  run  on  the  IBM 
System/360  machine  at  the  w .  R.  Church  Computer  Center, 
Naval  Postgraduate  School.  To  use  the  program,  a  deck  of 
plot  identification  cards  must  be  appended  following  the 
FORTRAN  Subroutine,  in  accordance  with  the  DSL/360  user's 
manual. 

The  labels  and  variable  names  used  in  the  program  are 
defined  in  the  Glossary  for  Frequency  Response  Programs, 
Appendix  D  of  this  thesis. 
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APPENDIX  B 


Program  fisting  for  Frequency  Response  of 
Ship  with  Passive  Tanks,  with 
FORTRAN  Subroutine  CALC 


The  Frequency  Response  program  was  written  in  Digital 
Simulation  Language  (DSL/360)  to  be  run  on  the  IB* 
Systea/360  machine  at  the  K.  R.  Church  Computer  Center, 
Naval  Postgraduate  School.  To  use  the  program,  a  deck  of 
plot  identification  cards  must  be  appended  following  the 
FORTRAN  Subroutine,  in  accordance  with  the  DSL/360  user's 
man  ual. 

The  labels  and  variable  names  used  in  the  program  are 
defined  in  the  Glossary  for  Frequency  Response  Programs, 
Appendix  D  of  this  thesis. 
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APPENDIX  C 


Program  Listing  tor  Frequency  ResDonse  of 
Activated  Ship/Tank  System,  with 
FORTRAN  Subroutine  CALC A 


The  Frequency  Response  program  was  written  in  Digital 
Simulation  Language  (DSL/360)  to  be  run  on  the  IBM 
System/360  machine  at  the  W.  R.  Church  Computer  Center, 
Naval  Postgraduate  School.  To  use  the  program,  a  deck  of 
plot  identification  cards  must  be  appended  following  the 
FORTRAN  Subroutine,  in  accordance  with  the  DSL/360  user's 
manual. 

The  labels  and  variable  names  used  in  the  program  are 
defined  in  the  Glossary  for  Frequency  Response  Programs, 
Appendix  D  of  this  thesis. 
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APPENDIX  D 


SYMBOL 
A0r A1 , A2 

BS 

BT 

CLOOP 
COM  P 
D 

G 

GIMAG 

GHEAL 

H 

IM 

IMH 

JS 

JST 

jr 

KS 

KT 

LOG  W 

MAG 

MAG  H 

MDB 

MDBS 

N 

OLOOP 

PHI 


Glossary  for  Frequency  Response  Programs 

MAIN 


MEANING 

A3  Coefficients  of  characteristic  equation 
Damping  coefficient  of  ship 
Overall  damping  coefficient  of  tank-fluid 
Closed-loop  transfer  function 
Compensator  (controller)  transfer  function 
Denominator  polynomial  of  open-loop  transfer 
function 

Open-loop  transfer  function 

Imaginary  part  of  open-loop  transfer  function 
Real  part  of  open-loop  transfer  function 
Closed-loop  transfer  function 

Imaginary  part  of  open-loop  transfer  function 
Imaginary  part  of  open-loop  transfer  function 
Roll  inertia  coefficient  of  ship 
Mutual-coupling  inertia  coefficient 
Equivalent  inertia  coefficient  of  tank-fluid 
Static  righting  coefficient  of  ship 
Equivalent  static  coefficient  of  tank-fluid 
Natural  logarithm  of  radian  frequency 
Magnitude  of  open-loop  transfer  function 
Magnitude  of  closed-loop  transfer  function 
Decibel  gain  of  open-loop  transfer  function 
Decibel  gain  of  closed-loop  transfer  function 
Numerator  polynomial  of  transfer  function 
Open-loop  transfer  function 
Phase  of  open-loop  transfer  function 
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Q 


Denominator  polynomial  of  open-loop  transfer 
function 

S  LaPlace  transform  variable 

SST  Ship-tank  response  in  feedback  path 


SUBROUTINES 


SYMBOL 

A 

AS 

BEAM 

BP 

BS 

BT 

BT1 

D 

DABS 

DCAL 

DISP 
El,  E2 


GM 

GRAV 

H 

11,12,13,14 

JS 


MEANING 

Cross-sectional  area  of  vertical  leg  of  U-tube 
Cross-sectional  area  of  horizontal  leg  of 
U-tube 

Transverse  dimension  of  ship  (beam) 

Velocity  head-loss  coefficient  of  pump 
Damping  coefficient  of  ship 
Overall  damping  coefficient  of  tank-fluid 
Skin-friction  damping  coefficient  of  tank 
fluid 

Perpendicular  distance  from  center  of  rotation 
of  ship  to  horizontal  cross-leg  of  U-tube 
Absolute  value  (magnitude)  of  DCAL 

Calculated  value  of  D  to  force  the 
mutual-coupling  inertia  coefficient  to  zero 
Displacement  of  ship 

Incremental  diffenence  between  current  and 
previously  calculated  values  of  ship 
parameters 
Metacentric  height 
Acceleration  of  gravity 
Height  of  fluid  column  in  vertical  leg  of 
U-tube 

Evaluation  of  integrals  as  described  in  main 
body  of  thesis 

Roll  inertia  coefficient  of  ship 
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JS  1 

jsr 

JT 

KR  AD 

KS 

KS1 

KT 

L 

LEN 
LfID  A 

n 

R 

RHO 

RHDG 

TJ 

TK 

TPHI 

WIDE 

WN2 

X 


Roll  inertia  coefficient  of  unmodified  ship 
Mutual-coupling  inertia  coefficient 
Equivalant  inertia  coefficient  of  tank-fluid 
Radius  of  gyration  of  the  mass  of  the  ship 
Static  righting  coefficient  of  ship 
Static  righting  coefficient  of  unmodified  ship 
Equivalent  static  coefficient  of  tank-fluid 
Longitudinal  dimension  of  a  vertical  leg  of 
the  D-tube 

Longitudinal  dimension  of  ship  (length) 

Design  ratio  X 

Number  of  iterations  of  design  procedure 

Lever  arm  of  U-tube 

Mass  density  of  tank-fluid 

Weight  density  of  tank-fluid 

Temporary  JS 

Temporary  KS 

Natural  period  of  roll  of  ship 
Transverse  dimension  of  vertical  leg  of  U-tube 
Natural  frequency  of  ship-tank  system,  squared 
Design  ratio  X 
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APPENDIX  E 


Waveslope  Spectrum  of  the  Ocean 

The  waveslope  spectrum  of  the  ocean  was  generated  from 
the  relationships  of  the  trochoidal  wave  model.  It  is 
intended  to  be  used  in  conjunction  with  the 
Pierson -Moskowitz  wave  power  spectra. 

This  appendix  includes  the  listing  of  the  DSL/360 
program  which  generated  the  spectrum,  and  listings  of  the 
spectra  for  sea  states  3,  h,  5,  and  6. 
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RESET  LABEL 

LABEL  WAVE SLOPE  SPECTRUM  FOR  H=i5.0  FEET 


HAVESL3PE  SPECTRUM  OF  THE  OCEAN 
SIGNIFICANT  HEIGHT  =  3.3  FEET 
********  (SEA  STATE  3)  ******** 


1.0000E-01 
1 . 0593E-0 1 

1.  1220E-01 

1.  1885E-01 
1. 2589E-01 

1.33  35E-0 1 
1. 4125E-01 
1.4962E-01 
1. 5849E-0 1 
1 .6788E-01 

1.  7783E-01 
1 . 8836E-01 
1. 9953E-01 

2.  1135E-01 

2.  2387E-01 
2. 3714E-0 1 
2. 5119 E- 01 
2.6607E-01 
2.8184E-01 
2. 9854E-01 

3.  16  23E-01 
3.3497E-01 

3.  54 8  IE-  0  1 
3. 7584E-01 
3.  98  1 1E- 0  1 
4.2170E-01 
4.4668E-0  1 
4.7315E-01 
5.0119E-01 
5. 3088E-Q 1 
5. 6234E-0  1 
5 . 9566E-0  1 

6.  30 96 E-  0  1 
6 . 6834E-0 1 

7.  0794E-0  1 
7.4989E-G1 
7. 9433E-01 
8.41 39E-0  1 

8.  9125E-01 
9.4406E-01 
1.0000E  00 
1.0593E 
1.  1220E 
1  .  1885E 
1. 2589E 
1 .3335E 
1.4125E 
1 .4962E 
1. 5849E 
1 .6788E 
1. 7783E 

1: 


mi 


00 

00 

00 

00 

00 

00 

09 

00 

00 

00 

n 


PHIMAX 

9.  3579E-03 
1 . 05  00E-02 
1.  17  8  IE-  02 
1.32 18E-02 
1 . 48  3  IE- 02 
1.6641 E-02 
1 . 867  IE-02 
2. 0950E-02 
2.35C6E-Q2 
2.63  74E-02 
2.95  9  2E- 02 
3.32  03  E-02 

3.  7254E-02 

4.  18  00  E-02 
4 . 69  00E- 02 
5.26  23  E-02 
5. 9044 E-02 
6.6249E-02 
7.4332E-02 

8. 34  02E-02 
9.3579E-02 
1.05  00  E-0  1 
1.  178  IE-01 
1.3218E-01 
1 . 48 3  IE- 0 1 
1.664 1E-0  1 
1. 8671E-01 

fcBJSfcSl 

2. 63 74 E-0 1 
2.9592E-01 
3.32  03 E-0  1 
3.7254E-01 
4 . 18  00  E-0  1 
4. 6900E-Q1 
5.26  23  E-0 1 

5.  9044E-01 
6.6248E-01 
7 . 4  3  32E- 0 1 

8 . 34  02 E-0  1 

9.35  78E- 0 1 
1.0500E  00 
1. 178  IE  00 
1.3218E  00 
1.4831E  00 
1 .6641E  00 
1.8671E  00 
2.0949E  00 
2.35C6E  00 
2.6374E  00 
2.9592E  00 

i-.mi  n 
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WAV2SL0PE  SPECTRUM  OF  THE  OCEAN 
SIGNIFICANT  HEIGHT  =  3.3  FEET 


******** (sea 

STATE  3)  ****** 

H 

PHIMAX 

2. 1135E 

00 

4 . 1 8  00  £ 

00 

2. 2387E 

00 

4. 690QE 

00 

2.3714E 

00 

5.2623E 

00 

2.5119E 

9° 

5. 9044E 

00 

2.6607E 

00 

6.6248E 

00 

2. 8184E 

00 

7.4331E 

00 

2. 9854 E 

00 

8.  1851  E 

00 

3.  16  23E 

00 

8.  18  5  1 E 

00 

3 . 34  96E 

00 

8.  185  1  E 

00 

3. 548  IE 

00 

8.  185  1 E 

00 

3.7584E 

00 

8. 1851E 

00 

3. 98  1 1 E 

00 

8. 185  1 E 

00 

4.2169E 

00 

8.  1 8  5 1  E 

00 

4. 4668E 

00 

8.  18  5  1 E 

00 

4.7315E 

00 

8. 1851E 

00 

5.0118E 

00 

8.  1851E 

00 

5.3088E 

00 

8. 1851E 

00 

5.6234E 

00 

8.  18  5  IE 

CO 

5.9566E 

00 

8.  1851E 

00 

6. 3095E 

00 

8. 185  IE 

00 
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WAVESL3PE  SPECTRUM  OF  THE  OCEAN 
SIGNIFICANT  HEIGHT  =  6.9  FEET 
******** ^ SEA  STATE  4) ******** 


1.  OOOOE-OI 
1. 0593E-01 
1.  1220E-0  1 
1. 1885E-01 
1. 2589E-01 
1 . 3335E-0 1 
1.4125E-01 
1.49  62E-0 1 
1. 5849E-0  1 
1 . 6788E-0 1 
1.  7783E-01 
1. 8836E-01 

1.  9953E-01 

2.  1135  E—  0 1 
2.  2387E-01 
2.37 14E-Q1 

2.  5119E-01 
2 . 66  07E-0 1 
2. 8184E-Q 1 
2. 9854E-01 

3.  16  23E-01 
3.3497E-01 
3 . 54  8 1 E-  0  1 


4.4668E-Q1 
4.7315E-01 
5. 0119  E- 01 
5. 3088E-0 1 
5. 6234E- 0  1 
5. 9566E-01 

6.  30  96E-01 
6 . 6834E-0 1 

7 .  07  94  E-  0  1 
7. 4989E-0 1 

7.  9433E-01 
8.41 39E-0 1 

8.  9125E-01 
9 . 44  Q6E-Q 1 
1.0000E  00 
1.0593E  00 
1.  1220E  00 

1:11111  88 
1.3335E  00 
1.4125E  00 

1:UH 1 88 

1.6788E  00 
1.7783E  00 

1:11181  QQ 


PHie 

I  AX 

1. 

95  66E- 

02 

2. 

19  54  E- 

02 

2. 

4633E- 

02 

2. 

7638E- 

02 

3. 

10  1  IE- 

02 

3. 

47  95  E- 

02 

3. 

9040E- 

02 

4. 

3804E- 

02 

4. 

9149E- 

02 

5. 

5146E- 

02 

6. 

1875E- 

02 

6. 

94  24  E- 

02 

7. 

76S5E- 

02 

8. 

74  00E- 

02 

9. 

8065E- 

02 

1. 

1003E- 

01 

1. 

2346E- 

01 

1. 

38  52  E- 

01 

1. 

5542E- 

01 

1. 

74  39  E- 

01 

1. 

9566E- 

Q  1 

2. 

19  54  E- 

01 

2. 

4633E- 

01 

2. 

7638E- 

01 

3. 

1011E- 

01 

3. 

4795E- 

01 

3. 

9040E- 

01 

4. 

3804  E- 

01 

4. 

9149E- 

01 

5. 

5146E- 

01 

6. 

1874E- 

Cl 

6. 

94  24  E“ 

01 

7. 

78S5E- 

01 

8. 

74  00  E- 

01 

9. 

8064E- 

01 

1. 

1003E 

00 

1. 

2346E 

00 

1. 

3852  E 

00 

1. 

1. 

1. 

5542E 
74  39  E 
9566E 

00 

88 

2. 

1954E 

00 

2. 

4633E 

00 

1: 

76  38  E 
10  1 1 E 

88 

3. 

47  94  E 

00 

3. 

9040E 

00 

4. 

4. 

9T48E 

88 

5. 

5 1  45  E 

00 

6. 

1874E 

CO 

1: 

94  24  E 
78S5E 

00 

CO 
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WAVESLO PE  SPECTRUM  OF  THE  OCEAN 
SIGNIFICANT  HEIGHT  =  6.9  FEET 
******** (SEA  STATE  4)  ******** 


2.1135E  00 
2.2387E  00 
2.3714E  00 
2.5119E  00 
2.6607E  00 
2.8184E  00 
2. 9854 E  00 
3.  1623E  00 
3.3496E  00 
3.  54 8 IE  00 
3.7584E  00 
3. 98 1 1E  00 
4.2169E  00 
4.4668E  00 
4.7315E  00 
5.0118E  00 
5.3088E  00 
5. 623UE  00 
5.9566E  00 
6.3095E  00 


PHIMAX 


8.1851E  00 
8.  185  IE  00 
8.1851E  00 
8.  185  IE  00 
8.  185 1 E  00 
8.1851E  00 
8.1851E  00 
8.1851E  00 
8.1851E  00 
8. 18 5  1 E  CO 
8.  IS 5 1 E  00 
8.  18 5 1 E  00 
8. 18 5 1 E  00 
8.1851E  00 
8.  18 5  1  E  00 
8. 185  IE  00 
8.1851E  00 
8. 18  5 1 E  00 
8. 1851 E  00 
8.185U  CO 
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WAVESLOPE  SPECTRUM  OF  THE  OCEAN 
SIGNIFICANT  HEIGHT  =  10.  FEET 
****** •♦(SEA  STATE  5) ******** 


1.  OOOOE-O  1 
1.0593E-01 

1.  122QE-0  1 

1.  1885E-01 
1. 2589E-0  1 
1 .3335E-01 
1.41 25E-  0  1 
1 . 4962E-0 1 
1. 5849E- 0  1 
1 .6788E-01 
1. 7783E-01 
1 . 8836E-01 
1. 9953E-0  1 

2.  1135E-01 

2.  2387E-Q1 
2 . 37 14 E- 0 1 

2.  5119E-01 
2.6607E-01 

2.  8184E-01 
2.9854E-01 

3.  1623E-01 
3.3497E-01 

3.  548  IE- 0  1 
3. 7584E-01 

3.  98 1 1E- 0  1 
4 . 21 70E-0 1 
4. 4668E-0  1 
4.7315E-01 

5.  01 19E-0  1 
5. 3088E-01 

5.  6234E-01 
5. 9566E-Q 1 

6.  30  96E-01 
6 . 6834E- 0  1 

7.  0794E-0  1 
7.4989E-01 
7.  9433E-01 
8.41 39E-0 1 
8.9125E-01 
9. 44Q6E-Q1 
1.000QE  00 
1.0593E  00 

1.  1220E  00 
1.1885E  00 
1.2589E  00 
1.3335E  00 
1.4125E  00 
1.4962E  00 
1.5849E  00 
1.6788E  00 
1.7783E  00 
1.8836E  00 
1.9953E  00 


PHIHAX 

2.  83 57E-02 

3.  18  17E-02 

3.  5700E-02 
4.0056E-02 
4.4943E-02 
5.0427E-02 
5.6580E-02 

6.  34  84  £-02 
7. 1230E-02 
7.9921E-02 
8. 96732-02 
1. 0062E-0 1 
1.  1289E-01 
1 . 26  67  E-0  1 
1.4212E-01 
1.5946  E-0  1 
1.7892E-01 
2 . 00  75  E-0 1 
2.25  25E-01 

3.  18  17E-01 
3.5699E-01 
4. 0056E-0 1 

4.  4943E-01 
5.04  27  E-0 1 
5.6580E-01 
6.3484E-01 

7.  1230E-01 
7. 992  IE-01 
8. 96  73 E-  0 1 
1.0061E  00 
1.128SE 
1.2667E 
1.4212E 
1. 5946E 
1.7892E 
2. 00 75 E 
2.2525E 


3.  18  17E 

3.  5699E 
4.0055E 

4.  4943E 
5 . 04  27  E 
5.6580E 
6.3483E 

7.  12  30E 
7.9921E 

8.  1851E 
8.  1851  E 
8.  18  5  1 E 


00 

00 

CO 

00 

00 

00 

oc 


8°0 

00 

CO 

00 

CO 

00 

CO 

00 

00 

00 

CO 

00 

CO 
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WAVESLOPE  SPECTRUM  OF  THE  OCEAN 
SIGNIFICANT  HEIGHT  =  10.  FEET 


******** (SEA 

STATE  5)  ****** 

W 

PHIHAX 

2 . 1 1 35E 

00 

8. 1851E 

00 

2. 2387E 

00 

8. 1851E 

00 

2.3714E 

00 

8. 1851  E 

00 

2.5119E 

00 

8. 185  IE 

00 

2.6607E 

00 

8. 1851  E 

00 

2. 8184E 

00 

8. 1851E 

CO 

2. 9854E 

00 

8.1851E 

00 

3.  1623E 

00 

8.  18  5  IE 

00 

3.3496E 

00 

8. 18  5  1  E 

00 

3. 5481E 

00 

8.  18  5  1 E 

00 

3.7584E 

00 

8.  18  5  1  E 

00 

3. 981  IE 

00 

8. 185  IE 

00 

4.2169E 

00 

8. 18  5  1  E 

00 

4.4668E 

00 

8. 185  IE 

00 

4.7315E 

00 

8. 1851E 

00 

5. 01 18E 

00 

8.  18  5 1 E 

00 

5.3088E 

00 

8.  1851  E 

00 

5. 6234E 

00 

8. 18  £  1 E 

CO 

5.9566E 

00 

8.  18  5 1  E 

00 

6.3095E 

00 

8. 185  IE 

00 
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WAVESL3PE  SPECTBOM  OF  THE  OCEAN 
SIGNIFICANT  HEIGHT  =  15.  FEET 
********( SEA  STATE  6) ******** 
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HAVESLDPE  SPECTRUM  OF  TEE  OCEAN 
SIGNIFICANT  HEIGHT  =  15.  FEET 
******** (SEA  STATE  6) ******** 


2.  1135E 
2. 2387E 
2.3714E 
2.5119E 
2.66  07S 
2. 8184E 
2.9854E 

3.  1623E 
3.34  96E 
3. 5481E 
3. 7584E 
3. 98  1 1 E 
4.2169E 
4. 4668E 
4.7315E 
5. 0118E 
5.3088E 
5.6234E 
5.9566E 
6.3095E 


00 

00 

00 

00 

00 

00 

00 

88 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 


PHIMAX 


8.  18  5 1 E 
8.  18  5  IE 
8.  185  1 E 
8.  18  5  IE 
8.  185 1  E 
8.  18 5  IE 
8 . 18  5  1  E 
8.  185  IE 
8.  1851E 
8.  185  1 E 
8.  18  5 1  E 
8.  18  5  1 E 
8. 1851E 
8.  18  5 1 E 
8. 1851  E 
8.  185  IE 
8. 18  5 1  E 
8.  1851E 
8.  18  5 1  E 
8.  1851E 


00 

CO 

00 

00 

00 

00 

00 

90 

00 

CQ 

00 

00 

00 

00 

00 

00 

00 

CO 

00 

00 
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